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We present the first results from a detailed analysis of photometric and spec- 
trophotometric data on the narrow-line Seyfert 1 galaxy Mrk 335, collected over 
a 120-day span in the fall of 2010. From these data we measure the lag in the 
HellA4686 broad emission line relative to the optical continuum to be 2.7 ± 
0.6 days and the lag in the H/3A4861 broad emission line to be 13.9 ± 0.9 days. 
Combined with the line width, the Hell lag yields a black hole mass, Mbh = 
(2.6 ± 0.8) X 10'' Mq. This measurement is consistent with measurements made 
using the H/3A4861 line, suggesting that the Hell emission originates in the same 
structure as H/3, but at a much smaller radius. This constitutes the first ro- 
bust lag measurement for a high-ionization line in a narrow-line Seyfert 1 galaxy 
and supports a scenario in which the He ll emission originates from gas in virial 
motion rather than outflow. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert 



INTRODUCTION 



Narrow-line Seyfert 1 galaxies (NLSls) are a subset of active galactic nuclei (AGNs 
that show narrower broad emission-line compo nents than typical Type 1 AGNs, as well 



as a number of o ther distinguishing properties ( lOsterbrock fc Poggd Il985l . iGoodrichI Il989 



Boiler et al.lll996l ). Explanations for their unique characteristics include the possibility tha t 
they are either low-inclination or high-Eddington rate accreters (or both - see lBorosonll201 ll ) . 
Substantial blue enhancements in high ionization lines such as C lvA15 49 and He llA1640 are 



apparently typical o f, although not restricted to, NLSl galaxies (e.g., [Richards et al. 



Sulentic et al. 



2011 



2002 



2000[). This m ay be evidence for material in a disk wind ( iRichards et al.l 



Leighly &: Moord 120041 ). If the blue enhancement is due to a wind, the use of high 



ionization emission lines to measure virial black hole masses (Mbh) in these objects may be 
problematic, as the method relies on the assumption that the emitting gas is in virial motion 
around the black hole. 



To investigate th e structure in this region, we turn to reverberation mapping ( iBlandford fc McKee 



19821 . iPetersoru 119931 ). This method has been extensively used to measure the physical size 
of the broad line region (BLR) in Type 1 AGNs, including NLSls. Reverberation mapping 
relies on the correlation between variations of the AGN continuum emission and the subse- 
quent response of the broad emission lines that are seen in Type 1 AGNs. By monitoring 
AGN spectra over a period of time, one can measure the radius of the emitting gas from the 
central source by observing a time delay between variations in the continuum and emission 
line fluxes. Assuming the gas is in virial motion, this radius can be combined with BLR gas 
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velocity dispersion estimates to obtain a measurement of Mrh- To date, this ra e thod has been 
applied to measure BL R radii in nearly 50 AGNs (e.g., 
2010l . lBentz et al.lboioh . 



Peterson et al.ll2004j . iDenney et al. 



If broad high-ionization lines like He ll A4686 are emitted from virialized gas near the 
black hole, we expect much shorter reverberation time lags for HellA4686 than for low 
ionization lines like H/3A4861 because both ionizatio n stratification a n d the line width require 
this gas to be much closer to the central source. iPeterson et al.l (120001 ) investigated the 
NLSl NGC4051 and detected very broad, blue-enhanced HellA4686 emission in the RMS 
spec trum. Unfortunately , their time resolution was inadequate to reliably measure a Hell 
lag. iDenney et al.l ( 120091 ) measured an improved H/3 lag for NGC 4051 using d a ta fro m their 
2007 campaign, but were similarly unable to recover a Hell lag. iBentz et al.l (|2010( ) report 
marginal detections of Hell lags in two NLSl galaxies observed in their 2008 observing 
campaign. These recent studies achieved similar high sampling rates, but the expected He ll 
lags for these targets are too short for a robust detection in these datasets. We have therefore 
been unable to test whether this emission originates in outflowing gas or is in virial motion. 

We included the NLSl galaxy Mrk 335 in a recent reverberation mapping campaign 
that will be described in detail elsewhere (Grier et al., in preparation). One goal of this 
high sampling rate program was to measure the reverberation lag for a high-ionization line, 
HellA4686, in this source. Here we present the Hell results, having obtained a high enough 
sampling rate to measure its short time delay, and show that the Mbh estimate from the 
high ionization He ll line agrees with those from low ionization lines. 



OBSERVATIONS 



In general, we follow the observational and analysis practices of IDenney et al.l ( 120101 ). 
which largely follows the analysis described by lPeterson et al.l (120041 ) . Details on observations 
and subsequent analysis techniques will be discussed in the accompanying work by Grier et 
al. 



2.1. Spectroscopy 

The majority of the spectra were obtained using the 1.3m McGraw-Hill telescope at 
MDM Observatory. We used the Boiler and Chivens CCD spectrograph to obtain 82 spectra 
over the course of 120 nights from 2010 Aug 31 to Dec 28. We used the 350 mm~^ grating 
to obtain a dispersion of 1.33 Apixel"^, with a central wavelength of 5150 A and overall 
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spectral coverage from roughly 4400 A to 5850 A. The sht was oriented north-south and set 
to a width of 5".0 and we used an extraction window of 12".0, which resulted in a spectral 
resolution of 7.9 A. Figure [T] shows the mean and root mean square (RMS) spectra of Mrk 
335 from the MDM spectra. 

We also obtained 7 spectra with the Nasmith spectrograph and SPEC- 10 CCD at the 
2.6m Shajn telescope at the Crimean Astrophysical Observatory (CrAO). We used a 3".0 slit 
at a position angle of 90°. Spectral coverage was from approximately 3900 A to 6100 A. 



2.2. Photometry 



We collected 25 epochs of F-band photometry from the 70-cm telescope at CrAO using 
the AP7p CCD at prime focus, covering a 15 ^ x15^ field of view. The flux was measured 
within an aperture of 15".0. See ISergeev et al.l ( l2005l ) for more details. 



We also obtained 19 epochs of IZ-band photometry at the Wise Observatory of Tel- Aviv 
University using the Centurion 18-inch telescope with a 3072 x 2048 STL-6303E CCD and 
a field of view of 75'x50'. For these data, we u sed the ISIS image sub t raction pack age rather 
than aperture p hotometry to measure flux es (lAlard fc Luptonlll998l : lAlardI |2000[ ) following 
the procedure of IShappee fc Stanekl ( 120111 ). 



3. LIGHT CURVES AND TIME SERIES ANALYSIS 

3.1. Light Curves 

The reduced sp ectra were fiux-calib rated assuming that the [O ill] A5007 emission line 



flux is constant (see iDenney et al.ll2010l for details on data processing). Emission-line light 
curves were created for both the MDM and CrAO data sets by fitting linear continua under- 
neath the H/3 and He ll lines and integrating the flux above them. H/3 fluxes were measured 
between 4910-5100 A in the observed frame, with the continuum interpolation defined by 
the regions 4895-4910 and 5215-5240 A. Hell A4686 fluxes were measured from 4660-4895 A, 
with the continuum defined from 4550-4575 and 4895-4910 A. The CrAO light curves were 
then scaled to the MDM light curve to account for the different amounts of [O ill] light that 
enters the slits due to differences in seeing, slit orientation, and aperture size. 

The continuum light curve was created by taking the average 5100 A continuum flux 
of the MDM spectra, measured between 5215-5240 A in the observed frame. This light 
curve was then scaled and merged with the other continuum and photometric light curves 
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with corrections for the host galaxy starhght in the different apertures (see iPeterson et al. 



19951 ) . The final continuum and emission-line light curves are shown in Figure [2j Light curve 



statistics are given in Table [TJ 



3.2. Time delay measurements 



For c omparison with pr e vious resu lts, we used the interpolat ion method originally de- 
scr ibed bv iGaskell &: Sparkd ( 119861 ) and lGaskell fc PetersonI ( 119871 ) which was later modified 
by IWhite fc PetersonI (Il994j ) and Peterson et al. (1998, 2004) to measure the time lag. We 
cross-correlated the two light curves with one another, calculating the value of the cross cor- 
relation coefficient r at each value of time lag. Figure |3] shows the resulting cross correlation 
functions (CCFs) for the light curves. U ncertainties in lags ar e calculated usi ng Monte Carlo 



simula tions that employ the methods of iPeterson et al.l ( 1l998l ) and refined by lPeterson et al. 



(120041 ). For each realization, we measure the location of the peak value of the cross correla- 
tion coefficient (Tpeak.ccp), and the centroid of the CCF (Tcent.ccp), calculated using points 
surrounding the peak. We adopt the mean rpeak,ccF and rcent,ccF from the Monte Carlo 
realizations for our delay measurements and the standard deviation as our formal uncertain- 
ties. Before the light curves were cross correlated, we removed the l ong-ter r a line ar upward 
trend that is clearly visible in all three light curves (see Figure IWelshI (119991 ) discusses 
the value in this practice of "detrending" the light curves, as the cross correlation function 
(CCF) tends to latch onto long-term trends unassociated with reverberation, often resulting 
in incorrect lags. From our cross correlation analysis, we measure Tcent,ccF(H/3) = 13.9 ± 
0.9 days and Tcent,ccF(Hell) = 2.7 ± 0.6 days. All lag measurements are listed in Table [2l 

Previous reverberation studies hav e relied o n these fairly simple cross correlation meth- 
ods to measure r. Recently, however, IZu et al.l ( 120111 ) discussed an alternative method of 
measuring reverberation time lags called the Stochastic Process Estimation for AGN Rever- 
beration (SPEAR) and demonstrated its ability to recover accurate time lags. The basic 
idea is to assume all emission-line light curves are scaled and shifted versions of the contin- 
uum light curve. One then fits the light curves using a dam ped random walk model (e.g. 
Kelly et al.ll2009l . iKozlowski et al.ll2010l iMacLeod et al.ll2010l ) and then aligns them to deter- 
mine the tim e lag. Uncerta inties in lags are computed using a Markov Chain Monte Carlo 



method (see IZu et al.l l201ll ). SPEAR is remarkably good at predicting time lags in data 



sets with relatively large gaps in the sampling. Using SPEAR, we successfully recover time 
lags for both the H/3 and He ll emission lines. We allowed SPEAR to automatically remove 
the linear trend and include any resulting uncertainties in the overall lag uncertainties. We 
measure rspEAR(H/3) = 14.0 ± 0.3 days and rspEAR(Hell) = 1.6;I^q;5 days, also reported in 
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Table |2j We see good agreement with the CCF results from the H/3 emission line, and while 
there is a small difference between SPEAR and CCF lags for the Hell line, they are still 
statistically consistent with one another. We suspect this small difference is due to the gap 
in data that is very close to the peak in the He ll light curve. We adopt the CCF values for 
our mass calcuations to allow comparison with previous reverberation efforts. 



3.3. Line width measurement and Mbh calculations 



Assuming that the motion of the H/3-emitting gas is dominated by gravity, the relation 
between Mbh, hue width, and time delay is 



M, 



BH 



fcTAV^ 
G 



(1) 



where r is the measured emission-line time delay, AV is the velocity dispersion of the BLR, 
and / is a dimensionless factor that accounts for the structure within the BLR. The BLR 
velocity dispersion can be estimated using the line width of the measured broad emission line 
in question. This width is usually characterized by either the FWHM or the line dispersion, 
(Tiine- We use a]\„p be cause there is evidence that it produces less biased Mbh measure- 
ments (jPetersonll201ll ). We measure the line width in the RMS spectrum, which eliminates 
contributions from the contaminating narrow components. We adopt an average value of 

< / >= 5.5 based on the assurnption that AGNs follow the s ame Mbh -cr* rel ationship as 
quiescent galaxies (jOnken et al.ll2004l ). This is consistent with IWoo et al.l (120101 ) and allows 
easy comparison wit h previous results, bu t is about a factor of two larger than the value of 

< / > computed by iGraham et al.l (120111 ) . 



To determine the best value of aune, we use Monte Carlo simulations following lPeterson et al 



(12004 ). The resulting line widths are given in Table |2l Using our measured values of rcent,ccF 
for the average time lag and an^e from the RMS spectrum as AV, we compute Mbh using 
both the H/3 and Hell emission lines. We measure Mbh = (2.7 ± 0.3) x 10'' Mq using the H/3 
emission line and Mbh = (2.6 ± 0.6) xlO'' Mq using Hell. 



4. DISCUSSION 

As discussed above, several studies involving NLSl galaxies have found indications of 
outflows in high-ionization lines in the form of enhanced flux on the blue side of the emission 
lines. Inspection of Mrk 335 spectra from the HST archive shows this enhanced blueward 
flux is present in the ClvA1549 line as well, but the HellA1640 line is blended with Civ, so 
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we cannot see if it too exhibits this blue enhancement. In fact, the shape of the HellA4686 
emission hne in the RMS spectrum of Mrk 335 (Figure [1]) shows red and blue shoulders 
that could be a signature of disk structure. To search for possible outflow signatures in 
the HellA4686 emission line, we divided it into red and blue components and integrated 
each component separately, creating two He ll light curves. We then cross-correlated the red 
and blue light curves with one another to see if there is any time delay between the two 
components. Cross correlation analysis yields a centroid lag Tcent = 0.4 ± 0.8 days. This 
is consistent with zero and thus presents no evidence for bulk outflows in the HellA4686 
emission of Mrk 335. The consistency of the Mbh measurements made using the Hell lines 
with those from H/3 are also suggestive of virial motion rather than outflowing gas. 



Previous reverberatio n measurements of Mrk 335 wer e made bvlKasse 



3aum et a. 



fll997l) 



andlPeterson et al.l (Il998[) and subsequently reanalyzed by lPeterson et al.l (|2004| ) and lZu et al 
J201lh . IZu et all J201lh report a time delay of 15.3+2:2 days and dune ~ 920 km s"^ for H/3, 
but were unable to make a robust He ll measu rement, as their average time sampling was 
on the order of 10 days. iPeterson et al.l (120041 ) measure Mbh = (1-4 ± 0.4) x 10^ M© from 
the H/5 emission line. Our Mbh measurements deviate from theirs by almost a factor of two. 
We suspect the difference in Mbh is due to the difference in the line width measurements 
between the two campaigns. The uncertainties in line width measurements and in the / 
factor are the main sources of uncertainties in reverberation Mbh measurements - when 
the light curves are well -sampled, the lag rn e asurements themselves have been shown to be 
remarkably robust (e.g. I Watson et al.l |20 111 ). iPeterson et al.l (120041 ) (Table 6) find that the 
virial products computed for an object using data from different epochs often differ from 
one another by as much as a factor of two (e.g. NGC 5548) and that the typical fractional 
error in the virial products is about 33%. Given our uncertainties in the / factor and the 
limitations in trying to accurately describe the BLR velocity field with a single line-width 
characterization, we probably cannot actually do better than about a factor of two or three 
in individual Mbh measurements. 

Ground-based reverberation campaigns in the past have been limited to objects with 
H/3 time lags that are expected to be less than a month or two due to both the finite length 
of the campaigns (which typically last 50-100 days) and the fact that most objects are only 
observable from the ground for only about half of the year. Measurement of longer time 
lags would require extended campaigns, which are difficult to schedule. If, as our evidence 
suggests, the Hell emission line is in virial motion around the black hole, we can use this 
emission line to measure Mbh in objects at higher redshifts, as expected He ll lags in many 
of these high- luminosity objects are short enough to measure in one observing season. 
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5. SUMMARY 

We have presented the first robust He II A4686 reverberation lag measurement in a NLSl 
galaxy. We also measure the B./3 time lag in this galaxy and compute Mbh using both 
emission lines. The Mbh measurements from He ll and H/3 are consistent with one another, 
suggesting that the gas producing the Hell emission resides in the same structure as that 
producing H/3 emission. While other high-ionization lines such as CiV show evidence for 
outflows, we do not see this in He ll, possibly because the He ll-emitting gas does not arise 
cospatially with gas producing Civ emission. This has practical implications for future 
reverberation efforts, as the He ll emission may allow us to more efficiently measure Mbh in 
objects at high redshift. 
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Fig. 1. — The flux-calibrated mean (top) and RMS (bottom) spectra of Mrk335 in the 
observed frame {z — 0.02579). The flux density is in units of 10~^^ erg cm' -2 k-\ The 
dashed lines show the spectra before the [O ill] narrow emission lines were removed; the 
solid line shows the spectra after the subtraction. 
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Fig. 2. — Complete light curves for Mrk 335 from our observing campaign. The top panel 
shows the 5100A flux in units of 10^^^ erg s~^ cm~^A~^, the middle panel shows the flux 
in the HeIlA4686 region in units of 10~^^ erg s~^ cm~^, and the bottom panel shows the 
integrated H/3A4861 flux, also in units of 10~^^ erg s^^ cm~^. Open circles denote observations 
from MDM Observatory and asterisks represent spectra taken at CrAO. Closed squares 
show the photometric observations from CrAO, and closed triangles represent photometric 
observations from WISE Observatory. Vertical dashed lines have been placed at two obvious 
features in the continuum to aid the eye. The vertical lines have been shifted by the measured 
Hell and H/3 lag values (2.7 days and 13.9 days, respectively) to aid the eye in identifying the 
correct lag values for each emission line. Dotted lines show the trends that were subtracted 
before performing the cross correlation analysis. 
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Fig. 3. — CCFs for the light curves. Dashed hnes represent the autocorrelation function 
(ACF) of the continuum hght curve, and solid hnes show the CCFs for the emission lines. 
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Table 1. Light Curve Statistics 







Sampling 




Mean 






Time 




Interval (days) 


Mean 


Fractional 






Series 


N 


(T) 


T 

median 


Flux 


Error 


F 

^ var 


-^max 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


5100 A 


133 


1.0 


0.95 


7.57 ± 1.01 


0.013 


0.129 


1.53 ±0.04 


HeiiA4686 


89 


1.5 


1.00 


7.65 ± 1.00 


0.016 


0.130 


1.55 ±0.04 


H/?A1861 


89 


1.5 


1.00 


5.71 ±0.53 


0.015 


0.091 


1.38 ±0.03 



*Column (1) lists the spectral feature, and column (2) gives the number of points 
in the individual light curves. Columns (3) and (4) list the average and median time 
spacing between observations, respectively. Column (5) gives the mean flux of the 
feature in the observed frame , and column (6) shows the mean fractional error that 
is computed based on observations that are closely spaced in time. Column (7) gives 
the excess variance, defined by 

F (2) 

where is the fiux variance of the observations, S'^ is the mean square uncertainty, 
and (/) is the mean observed fiux. Column (8) is the ratio of the maximum to 
minimum fiux in each light curve. 

*Continuum and emission-line fiuxes are given in 10~^^ erg cm~^A~^ and 10~^^ 
erg s"-*^ cm~^, respectively. 
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Table 2. B./3 and HellA4686 Time Series Results 



Parameter 
(1) 



T'cent.CCF^ 
Tpeak,CCF^ 
TSPEAR 

(Jiine (mean) 
FWHM (mean) 
anno (RMS) 
FWHM (RMS) 



H/3 
(2) 



13.9 ± 0.9 days 
13.8 ± 0.8 days 
14.0 ± 0.3 days 
1641 ± 12 km 
1363 ± 15 km s"^ 
1336 ± 51 km s"^ 
1149 ± 38 km s"^ 
(2.7±0.3)xlO^M( 



Hen 
(3) 



2.7 ± 0.6 days 
2.1 ± 1.2 days 

1.6 t'^ol days 
3465 ± 26 km s"^ 
3191 ± 571 km s^^ 
3001 ± 277 km s^^ 
7380 ± 1275 km s"^ 
(2.6±O.8)xlO^M0 



*A11 values are given in the rest frame of the object. 



